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Advances in knowledge of the structure of the photon and tests of perturbative 
QCD have been made using the increased luminosity with the ZEUS detector at 
HERA. Events with a low photon virtuality and two high transverse energy jets 
have been studied. Measurements of inclusive dijet and multijet production are 
herein compared to Next-to-Leading-Order calculations. In order to analyse the 
structure of the photon, dijet production of quasi-real and virtual photons and 
dijet production containing D* mesons were measured. 

1 Introduction 

The study of dijet photoproduction at HERA allows one to perform tests of perturbative 
QCD (pQCD) and to place constraints on the structure of the photon. To Leading Order 
(LO) two types of processes contribute to jet photoproduction [1] - direct and resolved 
photon processes (see Fig. 1). In direct photon processes, the photon interacts directly 
in the hard sub-process, whereas for resolved photon processes, the photon resolves into 
a source of partons, one of which participates in the hard sub-process. 








(a) (b) (c) 

Fig. 1: Examples of LO direct photon (a) and resolved photon (b) and (c) processes. 

The cross-section for resolved photon processes can be written in terms of the per- 
turbatively calculable 2 — > 2 scattering cross-section, da a b^cd-, as, 



da~, 



cd = 



fp/bi x pi Mp)/ 7 /a( a; 7' tfydo-ab^cd 



(1) 



ab p x i 

The proton's parton density function, f p /b, is experimentally constrained allowing the 
extraction of the parton density of the photon, / 7 / OI a currently poorly known quantity. 
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Comparisons of data with pQCD calculations will be covered first followed by anal- 
yses of the strucure of the photon. 



2 Tests of pQCD in Dijet Photoproduction 

The invariant mass of the dijet system, M JJ , is sensitive to the presence of new particles 
or resonances decaying into jets. The scattering angle in the centre-of-mass frame of 
the dijet system, cosO* , reflects the underlying parton dynamics, thereby providing a 
test of QCD. In direct photon processes, the dominant propagator is that of a quark in 
the s, t and u channels but is a gluon in the t channel for resolved photon processes. 
The presence of the different propagators in the two processes is reflected in the angular 
dependence of the cross-section; the spin — \ quark yields a (1 — IcosO*])^ 1 and the 
spin — 1 gluon yields a (1 — \cos8*\)~ 2 dependence in the cross-section. This predicted 
(by pQCD) behaviour of the cross-section dependence on direct and resolved processes 
has been observed at HERA for M JJ > 23 GeV [2]. Consequently any deviation from 
the pQCD predictions for the cross-section for higher dijet invariant masses would also 
be an indication of the presence of decays from new particles or resonances. 

Photoproduction events were defined by requiring Q 2 < 4 GeV 2 (corresponding to 
a median Q 2 « 10~ 3 GeV 2 ) and the photon-proton centre-of-mass energy in the range, 
134 < W < 277 GeV. Differential cross-sections as a function of M JJ and cosO* for 
dijet masses above 47 GeV and |cos6**| < 0.8 were measured. Jet finding was performed 
using the kr 1 clustering algorithm [3] requiring there to be at least two jets with 
> 14 GeV and —1 < rf et < 2. The measured cross-sections da / d\cos9* | and 
da/dM JJ using 41.3 pb _1 are shown in Fig. 2 and compared to pQCD calculations [4]. 
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Fig. 2: Differential cross-sections, da/d\cos6*\ (left) and da/dM JJ (right) compared to pQCD 
calculations. The ZEUS data points show statistical errors (thick bars) and statistical plus 
systematic errors (thin bars) with the error due to the energy scale displayed as a band. 



All measurements presented in this paper were performed using the fcy clustering algorithm 
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The NLO calculations agree well in shape with the measured distributions for both 
da / d\cos9* \ and da/dM J . The predictions using the GRV-HO photon structure func- 
tion are closer in magnitude to the data compared to those from GS96. 



3 Tests of pQCD in Multijet Photoproduction 



Multijet production provides a test of 
pQCD predictions beyond leading order and 
additionally tests extensions to fixed order 
theories such as parton shower models. The 
three jet system can be visualised in the 
centre-of-mass frame by considering Fig. 3. 
Of interest are the angles 8 3 and ^3. 8 3 
is the angle between the highest energy jet 
and the beam direction and is analogous to 
the scattering angle 8* from the previous 
section. ^3 is the angle between the plane 
containing the three jets and the plane con- 
taining the highest energy jet and the beam 
direction. 



THREE-BODY REST FRAME 

1 + 2 — 3 + 4 + 5 




Fig. 3: Illustration of the angles #3 and ip3 
for a particular three jet configuration. 
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By requiring Q 2 < 1 GeV 2 and a photon- 
proton energy range, 134 < W < 269 GeV, 
photoproduction events are selected. These 
events are then required to have at least two 
jets with E^ 1 > 6 GeV and a third with 
ESp > 5 GeV in a region of pseudorapid- 
ity, |r/-' et | < 2.4. These jet requirements in- 
troduce a bias in the angular distributions 
by excluding those jets produced close the 
beam-line. The criteria; M 3 j > 50 GeV, 
\cos9 3 \ < 0.8 and 2E 3 /M 3J < 0.95, reduce 
this bias. 

The three jet invariant mass cross-section, 
da/dM 3 j, is shown in Fig. 4 and compared 
with 0(aaf ) calculations from two pairs of 
authors, Harris & Owens [5] and Klasen & 
Kramer [6], showing good agreement with 
the data. The Monte Carlo models HER- 
WIG and PYTHIA describe the shape well 
but lie 20-40% below the data. In Fig. 5, 
the angular distributions, cos9 3 and ip 3 are 
shown. These distributions show a diffcr- 
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Fig. 4: Differential cross-section, da/dM^j, 
compared to pQCD calculations and Monte 
Carlo predictions. The ZEUS data points 
show statistical errors (inner bars) and sta- 
tistical plus systematic errors (outer bars) 
with the error due to the energy scale dis- 
played as a band. 

ence from those obtained from phase space 

displaying a sensitivity to the QCD matrix elements. The cos8 3 distribution is similar 
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to that of a Rutherford scattering form and is well predicted by the QCD calculations 
which account for the spin of the propagator. On consideration of the distribution in 
■03, one can see a tendency for the three jet plane to lie near the plane containing the 
beam and the highest energy jet as predicted by the coherence property of QCD. 
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Fig. 5: Distributions for cosOg (left) and V>3 (right). The statistical and systematic errors are 
as described previously. 



4 High Et Dijet Cross-sections in Photoproduction 

By choosing events with suitably high Et dijets, one can study the sensitivity of the 
cross-sections to the photon's parton distribution, the proton's being well constrained 
from previous experiments. This provides a measurement, in the range of the parton's 
momentum fractions, x 7 , with a higher scale to those in e + e~ data. The current pa- 
rameterisations of the parton density of the photon in the high x region have large 
differences due to uncertainties in experimental measurements. In previous dijet mea- 
surements [8], a large excess of the measured cross-sections for E T et > 6 GeV was seen 
over the NLO predictions for resolved enriched samples. This was attributed to pos- 
sible contributions from non-perturbative effects such as multiparton interactions [7]. 
This effect should be reduced for increasing E T et . To reconstruct x 7 experimentally, we 
define the observable: 



.r 



obs _ ^jets 



V E 31 

jets 



)< I 



e 



(2) 



7 2yE e 

where the sum runs over the two highest Ex jets and yE e is the initial photon energy. 
We can then define direct and resolved processes by cutting on x° bs , such that direct 
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enriched requires x° bs > 0.75. 

Cross-sections are measured in two 
kincmatical regions of the photon and 
proton centre-of-mass; the nominal re- 
gion, 134 < W < 277 GeV and the 
region 212 < W < 277 GeV, which en- 
hances the low i° bs events. The pho- 
ton virtuality was again required to be 
less than 1 GeV 2 . Asymmetric require- 
ments on the jet transverse energy [4] 

or /■:.;: ' ; 

in the region 
sen. 

Fig. 6 shows the uncorrected a;° bs dis- 
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Fig. 6: 

tion (dots) compared to HERWIG Monte Carlo 
predictions for direct only (shaded histogram) and 
direct plus resolved (open histogram). 
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tribution compared to HERWIG Monte 
Carlo predictions. There is good agree- 
ment in shape between the data and 

Monte Carlo without the requirement of multiparton interactions. For lower transverse 
energy jets (£J|f > 6 GeV [8]), there was a large excess at low x° hs which is not seen 
here. 

Evidence of the non-requirement of 
multiparton interactions to describe the 
data can also be seen the distribution 
of transverse energy flow around the 
jets. Consideration of Fig. 7 demon- 
strates good agreement between data 
and Monte Carlo predictions. For high 
values of 77, the region of disagreement 
in lower Et dijet studies [8], the agree- 
ment remains good. From this obser- 
vation, coupled to the measurement of 
x° hs , we conclude that no additional 
processes above the leading logarithm 
parton shower model is necessary to 
describe the event shapes. 

Differential cross-sections with re- 
spect to the pseudorapidity of the sec- 
ond jet in bins of the pseudorapidity of 
the first jet are here presented. For the 
entire region in W, da/drf^ is shown 
for the whole region in i° bs and for the 
direct enriched (x° hs > 0.75) region m 
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Fig. 7: Uncorrected transverse energy flow around 
the jets compared to HERWIG Monte Carlo pre- 
dictions for three pseudorapidity regions. The 
flows are measured using the energy deposits in 
a band of ±1 unit from the jet centre in <f> as a 



function of A77 = rf 



„jet 



Fig. 8. The data is compared to NLO calculations from Klasen et al. [6] and Harris 
et al. [9] for the GS96 photon structure function and the GRV-HO structure function 
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for Klasen et al.. Reasonable agreement between data and theory is seen in shape and 
magnitude for both regions of £° bs . However, differences of the order of the systematic 
errors are seen between the two structure functions. 
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Fig. 8: Measured da/drf 2 et in 134 < W < 277 GeV in three regions of r)[ et for all z° bs (circles) 
and for x° bB > 0.75 (triangles) compared to NLO calculations. The data has statistical errors 
(thick bars), the sum of statistical and systematic errors (thin bars) and a band due to the 
uncertainty in the energy scale. 
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Fig. 9: Measured dajdrf^ 1 in 212 < W < 277 GeV in three regions of r/{ e for all i* s compared 
to NLO calculations. The data has statistical and systematic errors as described previously. 



Photon Structure and Hadronic Final States at ZEUS 



7 



Fig. 9 shows da jdr^ 1 for the region of high W with all other cross-section defini- 
tons remaining the same. The cross-sections for the whole range of i° bs are shown 
and compared to the calculations from Klasen et al. with the the two previously men- 
tioned structure functions. From these we see that this high W region leads to an 
increased sensitivity to the photon structure function particularly in the 1 < rrl^ 1 < 2 
region. In the central region of pseudorapidity the measured cross-sections lie above 
the predictions. 

The uncertainty from possible differences between the jets from final state parti- 
cles and NLO partons due to higher order contributions and hadronizations is partly 
estimated using the leading logarithm parton shower Monte Carlo, HERWIG and 
PYTHIA. The estimator (da / 'drff ) parton/ '(da / 'dri J 2 et ) — 1 for the narrower W region, 
212< W < 277 GeV is shown in Fig. 10 where the partons are those after the parton 
shower process. The uncertainities are compared with the experimental uncertainitics 
for the measured da/drff . For low rf 2 et the estimator from the HERWIG model is 
large (smaller for PYTHIA) and comparable with the systematic uncertainties. The 
theoretical uncertainty is smaller than the experimental error in the other regions of 
rf et being similar for HERWIG and PYTHIA. 
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Fig. 10: Estimation of hadronization corrections, (da/dril et )parton/(da/drj^ et ) — 1. The parton 
level cross-section, (dcr/drf^ ) pa rton , is calculated using the partons after the parton shower 
in the leading logarithm Monte Carlos, HERWIG and PYTHIA. The errors are displayed as 
previously defined. 



Another uncertainty comes from the ambiguity of the renormalisation and factori- 
sation scale, n, and is estimated by Harris et al. [9] to be of the order of 10% by varying 
the scale; Et/2 to 2Et, for the cross-sections in Fig. 8. 
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5 Real and Virtual Photons in Dijet Production 



Whilst progress has recently been made in studying the parton distribution functions 
(PDFs) for quasi-real photons, little information exists for those of virtual photons with 
low statistics data from the PLUTO collaboration and studies of the photoproduction to 
deep inelastic transition region from the HI collaboration [10] being the only published 
results. The expectation for the HERA data is that the contribution to the cross- 
section from resolved photon compared to direct photon processes should decrease with 
increasing photon virtualities. Therefore measurements of the evolution of the resolved 
photon component with Q 2 are expected to constrain the virtual photon PDFs and test 
pQCD. 
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Fig. 11: Uncorrected x™ cas in three regions of photon virtuality compared to HERWIG Monte 
Carlo. Direct (shaded histogram) and the sum of direct and resolved (open histogram) are 
shown, where the ratio of the two was determined from a fit in each bin of Q 2 . 
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As an estimator of ir° bs , 



defined which is analagous to Eqn. 2 
but uses calorimeter quantities; yjs as 
an estimator for y, £jf 4 > 5 GeV 
and -1.125 < < e c ' as < 1.875. Fig. 11 
shows the £™ eas distributions compared 
to HERWIG Monte Carlo predictions 
for three bins of Q 2 . The normalistaion 
of direct and resolved processes were 
extracted from a two parameter fit to 
the measured distribution, being differ- 
ent for each range in Q 2 . The agree- 
ment in shape for the virtual photon 
data is good but is poor for the quasi- 
real photon data due to non perturba- 
tive effects as mentioned in the previous 
section. 

Using this data we can form the ra- 
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Fig. 12: The ratio, a(x° hs < 0.75)/a(x° hs > 0.75) 
as a function of photon virtuality, Q 2 . The points 
have statistical errors (inner bars) and the sum 
of statistical and systematic (outer bars) and the 
band represents the uncertainty in the energy 
scale of the calorimeter. 



tio of resolved enriched (x™ cas < 0.75) 
to direct enriched (x™ oas > 0.75) events and then compute the ratio of cross-sections, 
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cr(x° bs < 0.75)/cr(x° bs > 0.75), as a function of Q 2 . This ratio defined for dijet events, 
E J T et > 6.5 GeV and -1.125 < n < 1.875, with 0.2 < y < 0.55 is shown in Fig. 12. From 
this figure we can see a general decrease in the ratio with increasing Q 2 . 



6 Charm in Dijet Photoproduction 



Allied to the jets, a charm quark provides an additional hard scale (m c > Aqcd) 
yielding a more reliable perturbative calculation. Two approaches for the calculation of 
charm are currently available, the so-called massive and massless schemes. The massive 
approach assumes only the three lightest quarks to be active flavours in the proton and 
photon. In the massless approach, the charm quark is treated as an additional active 
flavour. Massless calculations therefore predict, for a given factorisation scale, a larger 
resolved component compared to massive calculations. A measurement of i° bs will 
provide a method for the possibilty of distinguishing the two schemes and probing the 
question of charm in the photon. ZEUS 1 996+97 

The charm quark itself cannot be di- 
rectly tagged, consequently events are 
selected containing reconstructed D*^ 
mesons with px{D*) > 3 GeV/c and 
\rj(D*)\ < 1.5 in dijet photoproduction 
events. Before measuring a cross-section 
in cc° bs , we first consider the energy flow 
about the jet axis, Fig. 13. The jet pro- 
files, calculated in the same way as de- 
scribed in Fig. 7, are split into low and 
high x° bs for three bins of rf et for recon- 
structed jets; E° ALjet > 4 GeV. The 
HERWIG Monte Carlo without multi- 
parton interactions provides a good de- 
scription of the data even in the forward 
region where there was a discrepancy for 
inclusive dijets of this energy [8]. One 
can also see that the direct photon only 
Monte Carlo cannot describe the data 
at low values of An as can be seen in 
the x° bs < 0.75, < 7] jet < 1 bin. This 
excess energy flow in the rear direction 
is consistent with there being a photon 
remnant. 
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Fig. 13: Uncorrected transverse energy flow with 
respect to the jet axis compared to HERWIG 
Monte Carlo (full histogram) and LO-direct pho- 
ton Monte Carlo only (dotted histogram). The 
distributions are separated into low and high 
x° bs for three bins in rf et . 

Fig. 14 shows da/dx° bs in the range Q 2 < 1 GeV 2 , 130 < W < 280 GeV for jets 
with 1 77 1 < 2.4, £^ et1 ' 2 > 7, 6 GeV and at least one D* meson, satisfying the criteria 
previously mentioned. 
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Fig. 14(a) shows the HERWIG Monte Carlo (normalised to the data) agreeing in 



shape with the measured cross-section. 
There is a peak at high x° bs consis- 
tent with LO-direct photon processes. 
There also exists a significant cross- 
section at low s° bs which cannot be 
described by LO-direct only and needs 
some component of LO-resolved pho- 
ton processes. At the given scale the 
LO-resolved component is dominated 
by charm excitation processes of the 
form, Fig. 1(c). The required LO- 
resolved contribution is 45±5 (stat.)% 
(compared to the HERWIG prediction 
of 37%) of which 93% is from charm 
excitation processes. 

Fig. 14(b) shows a comparison of 
the data with an NLO massive calcu- 
lation [11]. This does not describe the 
low x° bs measured cross-section whilst 
describing the high x° hs region. Upon 
choosing extremes of the scale fin and 
the charm mass m c , the prediction still 
fails to describe the data. 
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Fig. 14: Differential cross-section, da/dx° bB for 
dijets with an associated D* meson compared to 
HERWIG Monte Carlo predictions (a) and a mas- 
sive NLO calculation (b). The inner (outer) error 
bars represent the statistical (statistical plus sys- 
tematic) errors and the band represents the un- 
certainty in the energy scale. 



7 Summary and Conclusions 

Sufficient data has now being collected by ZEUS to allow the comparison of pQCD with 

precise measurements; these comparisons show good agreement for dijet and multijet 

dynamics. With the increased precision the data can now differentiate between different 

photon structure function parameterizations and can be used to place constraints on 

its form. „ „ 
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